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Long-slit medium-resolution spectra of the Galactic globular clusters (GCs) NGC 6229 
and NGG6779, obtained with the GARELEC spectrograph at the 1.93-m telescope of the 
Haute-Provence observatory, have been used to determine the age, helium abundance (Y), 
and metallicity [Fe/H] as well as the first estimate of the abundances of G, N, O, Mg, Ga, 

Ti, and Gr for these objects. We solved this task by comparing the observed spectra and the 
integrated synthetic spectra, calculated with the use of the stellar atmosphere models with 
the parameters preset for the stars from these clusters. The model mass estimates, TeS, and 
log g were derived by comparing the observed “color-magnitude” diagrams and the theoret¬ 
ical isochrones. The summing-up of the synthetic blanketed stellar spectra was conducted 
according to the Ghabrier mass function. To test the accuracy of the results, we estimated 
the chemical abundances, [Fe/H], log t, and Y for the NGC5904 and NGC6254 clusters, 
which, according to the literature, are considered to be the closest analogues of the two GCs 
of our study. Using the medium-resolution spectra from the library of Schiavon et ah, we 
obtained for these two clusters a satisfactory agreement with the reported estimates for all 
the parameters within the errors. We derived the following cluster parameters. NGC 6229: 
[Fe/H] = -1.65 dex, t = 12.6 Gyr, Y = 0.26, [a/Fe] = 0.28 dex; NGC 6779: [Fe/H] = 

-1.9 dex, t = 12.6 Gyr, Y = 0.23, [a/Fe] = 0.08 dex; NGC 5904: [Fe/H] = -1.6 dex, 
t = 12.6 Gyr, Y = 0.30, [a/Fe] = 0.35 dex; NGC 6254: [Fe/H] = -1.52 dex, t = 11.2 Gyr, 

Y = 0.30, [a/Fe] = 0.025 dex. The value [a/Fe] denotes the average of the Ca and Mg 
abundances. 
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of GCs are old (t = 7-15 Gyr). The clusters move 
with different velocities, at different distances rel¬ 


ative to the Galactic center, and have different 
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abundances. Objects with relatively high metal- 
licity occur closer to the Galactic center. They 
are called bulge clusters. The GCs with the low¬ 
est metallicity (—2.4 < [Fe/H] < —1.5 dex) oc¬ 
cur in the halo of our Galaxy. The most distant of 
these objects is NGC 2419. Its galactocentric dis¬ 
tance is equal to Rqc = 84 kpc which is 3.5 times 
greater than the distance Rqc = 24 kpc of the 
nearest dwarf spheroidal to the Milky Way: the 
Sagittarius dSph. The typical parameters of GG 
subsystems and their population characteristics 
are described, e.g., in Borkova and Marsakov [2]. 

The metallicities of the majority of GGs are 
derived from a comparison of stellar photome¬ 
try results and theoretical evolutionary tracks 
and isochrones. The brightest cluster stars are 
investigated with high resolution spectrographs. 
By no means all of the Galactic GGs have been 
studied spectroscopically. The clusters of our 
study (NGG6229 and NGG6779 (M59)) belong 
to such poorly studied objects. 

The basic properties of the studied objects are 
presented in Table 1 along with the parameters of 
the objects we chose from [1]; these two objects, 
NGC 5904 (M5) and NGC 6254 (MIO), are their 
close analogues in terms of metallicity, age, hor¬ 
izontal branch type, and other parameters. We 
took the spectra of these GCs from Schiavon et 
al. [3]. All four GCs under consideration are mas¬ 
sive, dynamically evolved objects with high stel¬ 
lar density in the center. 

NGC 6229 is one of the most distant GCs of 
the outer halo. The photometry was performed 


and the “color-magnitude” diagram (G-M) was 
plotted for this cluster by several authors during 
1986-1991 (see references in [4]). The first deep 
photometry of the cluster core and periphery, 
including the full horizontal branch (HB) and 
the asymptotic giant branch, was performed by 
Borissova et al. with the 2-m telescope at NAO 
Rozhen, Bulgaria [4-6]. In these papers and also 
in [7], the HB structure and the variable stars 
in NGC 6229 were investigated; the position of 
the main sequence (MS) turnoff point was also 
determined. Borissova et al. [5] were the first to 
note the similarity in the metallicity and age of 
NGC 6229 and NGC 5904. A common property 
of these two clusters is the presence of a consid¬ 
erable number of RRLyrae-type variable stars. 
The blue stars brighter than the MS turnoff point 
(blue stragglers) were studied in [8] based on the 
Hubble Space Telescope images. 

NGC 6779 is a low-metallicity globular clus¬ 
ter of the halo. It is comparatively poorly stud¬ 
ied because of its close proximity to the Galactic 
plane. Hatzidimitriou et al. [9] presented a deep 
ground-based photometry of this cluster, and 
Sarajedini et al. [10] and Dotter et al. [11] per¬ 
formed the stellar photometry based on the Hub¬ 
ble images. The authors estimated the metallic- 
ity and age of the cluster using the Dartmouth 
isochrones [12]. Unlike NGC 6229, NGC 6779 has 
a low metallicity (see Table 1) and an extended 
blue horizontal branch. 

NGC 5904 is the closest GC to us that is 
located far from the Galactic plane. That is 
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Table 1. The key parameters of the studied and comparison globular clusters according to [1]: heliocentric 
distance Dist©; Galactocentric distance Distcc, assuming Distcc = 8 kpc for the Sun; absolute V'-band magnitude 
in the Johnson-Cousins photometric system; central y-band surface brightness; color excess; tidal radius rr, core 
radius rc; central concentration c = log(rt/re); central luminosity density logarithm 


Parameter 

NGC 6229 

NGC 6779 

NGC 5904 

NGC 6254 

Dist©, kpc 

30.4 

10.1 

7.5 

4.4 

Distcc, kpc 

29.7 

6.2 

9.7 

4.6 

Mv, mag 

-8.05 

-7.38 

-8.81 

-7.48 

/iy, mag arcsec”^ 

16.99 

18.06 

16.05 

17.69 

E{B — V), mag 

0.01 

0.20 

0.03 

0.28 

n, pc 

47.6 

25.15 

61.96 

27.49 

Tc, pc 

1.15 

1.08 

0.91 

1.10 

c 

1.61 

1.37 

1.83 

1.40 

logLcent, Lq/pC^ 

3.41 

3.26 

3.91 

3.5 


why it is an excellent object for studying the 
abundances of individual GC stars, and for plot¬ 
ting deep C-M diagrams. Its exact photome¬ 
try serves for testing the stellar evolution models. 
The cluster includes a large number of variable 
stars. It has an enormous space velocity and a 
very eccentric orbit [13]. Two millisecond pul¬ 
sars have been found in it, one of which is as¬ 
sociated with a massive neutron star [14]. The 
nature of these objects is not clear yet. Coppola 
et al. [15] estimated the distance to the cluster 
by the “period-infrared magnitude” relation for 
the RRLyr-type stars: (M —m)o = 14.44±0.02. 
Similar distance estimates were obtained by ap¬ 
proximating the evolutionary sequences of the 
cluster with theoretical tracks using the Hubble 
photometry [16, 17]. 

NGC6254 is one of the nearby clusters of the 
Galactic halo for which deep Hubble photome¬ 


try is available. In particular, the photometry 
results allow one to determine the percentage of 
binaries as a function of distance to the center 
of NGC6254 [18]. This, in turn, made it possi¬ 
ble to explain the low level of mass segregation 
in the cluster without invoking the hypothesis 
of an intermediate-mass black hole in the clus¬ 
ter center [19]. The deep photometry performed 
on Hubble images and estimates of the distance 
and Galactic absorption for NGC6254 were also 
presented in [10]. 

The integrated spectra of the two nearby clus¬ 
ters, NGC5904 and NGC6254, were included 
in the spectral library of Schiavon et al. [3], 
which was used by different authors for numer¬ 
ous studies afterwards. The deep C-M di¬ 
agrams for the clusters NGG6779, NGG5904, 
and NGG6254 [10, 16] were widely used for the 
investigation of the cluster evolutionary status 
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and testing the stellar population models. Van- 
denBerg et al. [20] determined the age, [Fe/H], 
and helium abundance index for NGC6779, 
NGC 590, and NGC 6254 applying their Victoria- 
Regina models. 

Nowadays, with the development of observa¬ 
tional and computing techniques for the analysis 
of the individual stars in GCs and also for the de¬ 
termination of their parameters, the method of 
synthetic spectra, calculated using stellar model 
atmospheres, is often used. By varying the abun¬ 
dances of different chemical elements, the de¬ 
viation of the theoretical spectra from the ob¬ 
served ones is minimized. The used spectral 
range should be large enough to include both 
a sufficient number of lines of the same element 
with different degrees of ionization and numerous 
lines of various elements. To obtain the fullest 
and most accurate data, one should obviously 
use spectra with the highest possible resolution. 
However, obtaining such spectra with a suffi¬ 
ciently high S/N ratio requires much observing 
time with large telescopes, and is possible only 
for the brightest stars. The analysis of the inte¬ 
grated light of the clusters allows us to efficiently 
process the data obtained with 1-2-m telescopes 
for all the objects in the Galaxy, and the data 
from large telescopes—for extragalactic clusters. 
This approach is adopted in this paper for the 
four above-mentioned clusters. It allowed us to 
obtain valuable data for the comparison of the 
photometric and spectral stellar evolution mod¬ 
els which is the essential basis for understanding 


the properties of the stellar populations in other 
galaxies. 

In Section 2 the observations and their reduc¬ 
tion techniques are described. In Section 3 we 
present the method of modeling the integrated 
cluster spectra in accordance with the data on 
the isochrones and the luminosity function, and 
of determining the chemical abundances. In Sec¬ 
tion 4 the results of the investigation are dis¬ 
cussed, and the conclusions are formulated in 
Section 5. 

2. OBSERVATIONS AND REDUGTION 

We have observed the clusters NGC 6229 
and NGC 6779 with the 1.93-m telescope of the 
Haute-Provence observatory for a comparative 
analysis of their stellar populations and the char¬ 
acteristics derived from the study of the C-M 
diagrams. The observation log is presented in 
Table 2. Long-slit medium- resolution spectra 
(5.5' X 2") were obtained with the CARELEC [15] 
spectrograph. We used a 300 lines/mm grating, 
which provided a dispersion of about 1.78 A per 
pixel and a ~ 5 A(FWHM) spectral resolution in 
the working spectral range of 3700-6800 A. Cal¬ 
ibration exposures with HeNe lamps were taken 
at the beginning and end of each night. The ob¬ 
servations were carried out in adverse weather 
conditions with cirrus and intermittent clouds. 
The average seeing was 2.5-3.5". However, the 
lack of illumination from the moon allowed us to 
obtain the CCD images of required quality dur- 
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Table 2. Log of spectroscopic observations 


Object 

Date 

Exposure, 

s 

Slit 

position, deg. 

NGC6229 

09.07.10 

600 

0 



600 

90 

NGC 6779 

09.07.10 

300 

95 



300 

99 


11.07.10 

2 X 900 

0 


ing clear sky periods. Table 2 shows that the 
observations were carried out with several fixed 
slit positions. In all cases the slit was oriented 
so as to capture not only the central, brightest 
part of a GC but also the neighboring regions 
highly populated with stars. Thus, the result¬ 
ing 2D image permitted us to compare the en¬ 
ergy distributions and individual spectral lines 
of the central region, unresolved into individual 
objects, and the neighboring stars. As a result, 
we could safely eliminate the background objects 
that could distort the GC abundances from the 
one-dimensional sum spectrum. Note that other 
authors usually use a scanning technique with 
a moving slit to obtain the integrated GC spec¬ 
trum (see, e.g., [3]). This technique enhances the 
probability of background objects entering into 
the resulting spectrum, especially in the case of 
GCs close to the Galactic plane. 

Primary data reduction was conducted us¬ 
ing MIDAS [21] and IRAF^ software packages. 
First, we filtered out cosmic ray hits with the 
filter/cosmic code in MIDAS environment. 

^ http://iraf.noao.edu 


Hot pixels were masked. The standard procedure 
of spectra reduction was then conducted: bias 
substraction, flat-field correction. The dispersion 
relation provided an accuracy of wavelength cal¬ 
ibration of about 0.8 A. One-dimensional spec¬ 
tra were extracted with the IRAF apsum proce¬ 
dure. The final signal-to-noise ratio per pixel 
at the center of the spectral range of the ob¬ 
tained spectra was S/N ~ 130. The contin¬ 
uum in the one-dimensional spectra was ap¬ 
proximated with the MIDAS filter/maximum 
and f ilter/smooth (running median) functions. 
The analysis of the correspondence of the ob¬ 
served and theoretical spectra was conducted in 
the Origin 6.1 graphical environment. 

3. METHOD OF DETERMINATION OF 
THE CHEMICAL ABUNDANCES AND 
EVOLUTIONARY PARAMETERS OF 
GLOBULAR CLUSTERS 

The method used in this paper was first pre¬ 
sented by Sharina et al. [22]. The method applies 
not only to spectral but also to photometric ob¬ 
servational and theoretical data. It is based on 
a comparison of the observed GC spectra and 
model spectra derived by summing up the indi¬ 
vidual synthetic spectra of stars with different 
masses, log g, and Tgff. The stellar parameters 
are set by the isochrone best corresponding to the 
cluster C-M diagram. The method for estimat¬ 
ing the abundances and evolutionary parameters 
described in this section can be applied to any 
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globular cluster for which deep stellar photom¬ 
etry and a long-slit spectrum containing the in¬ 
formation on all the stars of the cluster are avail¬ 
able. The spectra of the dense central GC re¬ 
gions are best suited for such an analysis. Our 
method is applicable to spectra with a resolution 
oi R > 2500, a rather high signal-to-noise ratio 
{S/N > 100), and a spectral range of no less than 
1500 A. 

3.1. Determination of Masses, log g, and Teg 
for Cluster Stars 

The synthetic blanketed spectra of stellar at¬ 
mospheres are calculated according to the phys¬ 
ical stellar parameters derived from the compar¬ 
ison of the observed C-M diagrams of the clus¬ 
ters and the theoretical isochrones from Bertelli 
et al. [23].^ For our work we chose the isochrones 
obtained by the Padova team. This choice was 
determined by the fact that the authors included 
the evolutionary stages of the HB and the asymp¬ 
totic giant branch in the models and also used 
a wide range of parameters: metallicity [Fe/H], 
age t, and helium abundance Y. For the clus¬ 
ters NGC6229 and NGC6779 we used the stel¬ 
lar photometry data from Piotto et al. [24] and 
Sarajedini et al. [10]. 

When selecting a theoretical isochrone which 
best fits the cluster C-M diagram, we vary five 
parameters: age t, specific helium abundance Y, 


metallicity [Fe/H], color excess E{B — V), and 
distance to the cluster. As initial parameters, we 
used data from the following references: Harris’s 
compiled catalog [1] and Vandenberg et al. [17], 
who determined the parameters based on stellar 
photometry of Sarajedini et al. [10] and Victoria- 
Regina [20] models. The initial metallicities and 
ages for a number of clusters were adopted from 
Dotter et al. [11], whose results were derived with 
the use of isochrones from Dartmouth [12]. 

Isochrone selection is carried out in the fol¬ 
lowing way. 


(1) The position of the basic evolutionary se¬ 
quences: main branch, subgiant branch, red 
giant branch, and horizontal branch, is com¬ 
puted from the photometric data as a run¬ 
ning median with a step of 0.1’ along two co¬ 
ordinates (color and magnitude). Stars with 
magnitude measurement errors larger than 0.2’ 
should previously have been excluded from the 
photometry table. Those are the faint objects. 
Usually, the percentage of their detection on 
Hubble images does not exceed 50%. Gonse- 
quently, the fiducial sequence of the C-M dia¬ 
gram is plotted, which describes the position of 
the basic evolutionary sequences of a cluster. 


(2) Next, the theoretical isochrone minimally 
deviating from the fiducial sequence is se¬ 
lected. This is done by y-square minimization: 

N r 


= 


E 

i=l 


fo I'm 


^lo 


^ http://stev.oapd.inaf.it/YZVAR/ 


where/o, crj^ are the magnitudes of the fiducial 
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Figure 1. The isochrone log t = 10.1, Y = 0.26, Z — 0.0004 for NGC6229. The circles of different sizes 
denote the points (in the circle centers) which contribute less than 1%, from 1% to 3%, and more than 3% to 

the total flux. 


and their errors, and Im are the corresponding 
isochrone values. Note that we do not inter¬ 
pret theoretical isochrones but merely use the 
available set [23]. 

Variations of the Galactic absorption in the 
direction of a GC and the distance to a cluster do 
not influence the inclination, form, and relative 
position of different evolutionary sequences, but 
shift the C~M diagram as a whole. Variations of 
[Fe/H], f, or V, on the other hand, influence the 
position and form of the isochrones (see, e.g., the 
explanation in [23]). The evaluation of the basic 


evolutionary parameters from the C-M diagram 
is conducted taking into account these variations. 
In the following six paragraphs, we explain the 
principles of determining [Fe/H], t, and Y. 

Selection of [Fe/H] (t, Y are constant) is 
performed by the inclination of the red giant 
branch [25] as the basic criterion. The inclina¬ 
tion of the branch increases with the increase of 
[Fe/H], and the luminosity of the stars at the 
top of the red giant branch falls more and more 
rapidly at [Fe/H] > —0.8 dex. Moreover, with 
increasing metallicity, the whole isochrone shifts 




toward the red region, the luminosity of the stars 
in it decreases, and fewer hot blue stars appear 
on the HB. 

Selection of the age value ([Fe/H] and Y are 
constant) is done using the temperature and lu¬ 
minosity of the MS turnoff stars. With increasing 
age, the MS turnoff point position shifts toward 
cooler and fainter stars and is additionally de¬ 
pendent on the GC chemical abundance. The 
luminosity of HB stars is mainly determined by 
metallicity and helium abundance. There are two 
basic and common empirical approaches to rela¬ 
tive age estimation for GCs: a vertical one [26- 
28], which takes into account the difference in 
the luminosities of the MS turnoff point and the 
zero age HB {6V = Vro “ Fzahb)) and a hori¬ 
zontal one [29], based on the difference in tem¬ 
perature (color) between the MS turnoff point 
and the base of the red giant branch 5{B — V) = 
{B — F)rgb — {B — H)to- 

In [27], an approximate ratio between dV, the 
abundances of helium Y and heavy elements Z 
on the one hand, and the logarithm of age on 
the other hand is given. As a consequence of the 
strong dependence of the HB star luminosities on 
Y, and the MS turnoff point on Z, the vertical 
method provides only the relative ages for GGs 
of similar metallicity and Y. For the estimation 
of absolute ages, the exact values of the helium 
mass fraction are needed. 

The essence of the horizontal method is in 
measuring the color difference between the MS 
turnoff point, sensitive to age, and the red giant 


branch, sensitive to metallicity. This method is 
also relative. The compared clusters should have 
the same metallicities and [a/Fe]. 

It is necessary to note that both the verti¬ 
cal and horizontal methods do not require knowl¬ 
edge of absolute distances to the observed objects 
and of the light absorption in their direction. 
There are numerous modifications of these two 
approaches (see, e.g., [30, 31]. The results of age 
estimation for the same cluster using different 
methods vary greatly. Modern researchers use a 
maximum likelihood method for isochrone selec¬ 
tion for the C-M diagram of a cluster as a whole. 
We use a similar method. Estimation of the MS 
turnoff point magnitude is rather difficult and 
often leads to an age inaccuracy of about 10%. 
That is why the approximation by an isochrone of 
the C-M diagram as a whole with all its bends 
is much more precise. As we use not only the 
photometric data characterizing the stellar lumi¬ 
nosities and temperatures but also spectroscopic 
data on the detailed abundances, further devel¬ 
opment of our method could allow us to estimate 
the absolute ages of the clusters. 

Variation ofY influences mainly the luminos¬ 
ity and temperature of HB stars. Old GGs have 
many blue stars on the HB and often a sequence 
of low-luminosity stars with temperatures higher 
than those of the HB stars [32]. 

After the selection of [Fe/H], t, and Y using 
the C~M diagram, the spectrum with model pa¬ 
rameters corresponding to the selected isochrone 
is calculated. Further refinement of the val- 
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ues [Fe/H], t, and Y is obtained with the use 
of spectroscopic data (see detailed explanation 
in paragraph 3.2). The long slit is adjusted to 
the cluster’s center. Consequently, the spectrum 
of compact GCs includes the information on all 
their stars. Provided that there is no contribu¬ 
tion from bright background stars, the integrated 
spectrum allows us to match the data derived 
from the C-M diagram with the results of syn¬ 
thetic spectra obtained with stellar atmosphere 
models. For calculating the spectra, the num¬ 
ber of isochrone points was optimized by elimi¬ 
nating the points contributing less than 0.5% to 
the total luminosity of a cluster. To estimate 
the contribution of individual points to the to¬ 
tal light of a cluster, the flux in the continuum 
at the wavelength A = 5000 A, calculated for 
the model atmosphere with the parameters of a 
given point (Tefj, log g and [Fe/H]), is multiplied 
by the squared radius of the star, the weight of 
the point in the full mass interval, and the cur¬ 
rent value of the mass function. The flux derived 
this way is further divided by the total flux from 
all the points of an isochrone. Figure 1 shows 
isochrone points which contribute to the total 
light of NGC6229 by less than 1%, from 1% to 
3% and more than 3% in circles of different sizes. 

The derived results of isochrone selection us¬ 
ing the C-M diagram and the spectra for the 
globular clusters under study are shown in the 
left panels of Figs. 2-5. On the right-hand 
side, isochrone positions with varied parame¬ 
ters are presented. For NGC6229 the varia¬ 


tion of the isochrone position with increasing 
Y is shown, for NGC6779—with increasing age, 
and for NGG 5904 and NGC 6254—with the vari¬ 
ation of Y and age. For the last two clus¬ 
ters, isochrones with parameters close to those 
estimated in [29] are presented. Distance and 
E(B — V) were changed according to variations 
of the above mentioned parameters. All GGs 
have certain stars and sequences which cannot 
be described by isochrones. They are hot blue 
HB stars, blue stragglers, certain bright stars de¬ 
viating from the red giant branch (background 
objects or asymptotic giant branch stars). In 
Section 3.3.2 we discuss the significance of the 
listed deviations during the analysis of the corre¬ 
spondence of theoretical spectra calculated with 
the specified isochrones to the observed ones. 
Figs. 6-9 show the spectral regions with hy¬ 
drogen lines calculated with the parameters of 
isochrones which are used in Figs. 2-5. As the 
hnal values of [Fe/H], t and Y, we adopted those 
that best describe the observed spectra. 

The results of isochrone selection for the thor¬ 
oughly studied clusters NGC 5904 and NGC 6254 
(Figs. 4 and 5) agree with the quantities from 
the literature (Table 3) [1, 17]. The differences 
in age estimates do not exceed 1 Gyr. However, 
the [Fe/H] value is regularly underestimated by 
0.1-0.3 dex for all the clusters except NGC 6779. 
For the two tested clusters, the helium abun¬ 
dance was found to be higher than in the pho¬ 
tometric results of VandenBerg et al. [17]. The 
a-element abundance for NGC 6254 that we mea- 
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Figure 2. Comparison of the C-M diagram 
V~{B — V) for NGC6229 [24] with the isochrones 
from Bertelli et al. [23]. The left panel shows the 
isochrone log t = 10.1, F = 0.26, Z = 0.0004 used 
in the GC spectrum modelling. On the right, the 
isochrone log t = 10.1, Y = 0.3, Z = 0.0004 (see 
Section 3.1) is shown for comparison. 


Figure 3. Comparison of the C-M diagram 
I-{y — I) for NGC6779 [10[ with the isochrones 
from [23[. The left panel shows the isochrone 
log t = 10.1, Y = 0.23, Z — 0.0001 used in the 
GC spectrum modelling. On the right, the 
isochrone log t = 10.15, Y = 0.23, Z = 0.0001 is 
shown for comparison. 



V-l 


V-l 



V-l V-l 


Figure 4. Comparison of the C-M diagram 
I-{y — I) for NGC5904 [24[ with the isochrones 
from [23[. The left panel shows the isochrone 
log t — 10.1, Y = 0.30, Z = 0.001 used in the GC 
spectrum modelling. On the right, the isochrone 
log t = 10.15, Y = 0.26, Z = 0.001 (see Section 
3.1) is shown for comparison. 


Figure 5. Comparison of the C-M diagram 
I-{y — I) for NGC6254 [10[ with the isochrones 
from [23[. The left panel shows the isochrone 
log t = 10.05, Y = 0.30, Z = 0.0004, used in the 
GC spectrum modelling. On the right, the 
isochrone log t = 10.15, Y = 0.26, Z = 0.0004 is 
shown for comparison. 
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Figure 6. The fit of the isochrone parameters to the observed spectrum of NGC6229 (the profiles of the 
H/3, H 7 , H(5 lines). The three upper panels show the applied model log t = 10.1, Y = 0.26, Z = 0.0004 (see 
the corresponding C-M diagram on the left of Fig. 2). The three lower panels show the model log t — 10.1, 
Y = 0.30, Z = 0.0004 (the right-hand side of Fig. 2) for comparison. 


sured turned out to be 0.2-0.3 dex lower than the 
value reported in the literature. On the whole, 
our errors of the determination of the main clus¬ 
ter parameters do not exceed 20%. 

However, we note the common problem (in 
the literature and ours) of describing the clus¬ 
ter C-M diagram using the evolutionary tracks 
and isochrones. Actually, none of the presented 
isochrones describes the C-M diagram of the 
cluster in all details. We selected the isochrones 
that provide the best agreement between the ob¬ 
served and theoretical spectra calculated using 
the distribution of stars by mass, log g, and Tgfj 


corresponding to these isochrones. Note that be¬ 
sides the above-mentioned particular stars with 
noticeable deviations from the isochrone, in some 
other cases (see, e.g.. Fig. 4), the spectrum anal¬ 
ysis implies a small shift of the MS turnoff point 
toward higher temperatures or a change of the 
HB position. For NGC5904 such changes may 
be caused by the increased percentage of hot 
blue HB stars, the EHB (Extended Horizontal 
Branch), and/or blue stragglers caught in the 
spectrograph slit. 
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Figure 7. The fit of the isochrone parameters to the observed spectrum of NGC6779 (the profiles of the 
H/3, H 7 , H(5 lines). The three upper panels show the applied model log t = 10.1, Y = 0.23, Z = 0.0001 (see 
the corresponding C-M diagram on the left of Fig. 3). The three lower panels show the model log t = 10.15, 
Y — 0.23, Z — 0.0001 (the right-hand side of Fig. 3) for comparison. 


3.2. Modelling of Stellar Spectra 

The number of stars with a given mass is 
calculated according to the Chabrier mass func¬ 
tion [33]. For obtaining the synthetic spectra 
of plane-parallel, hydrostatic stellar atmospheres 
for the derived isochronic parameters (Tefj, log^f, 
[Fe/H]), we use the CLUSTER [22] software. 
Model atmospheres are computed by interpo¬ 
lating the grids of Castelli and Kurucz [34] ac¬ 
cording to the technique described in [35]. For 
normalization we model two spectra simultane¬ 
ously: (1) with the allowance for atomic lines and 


molecular bands in the spectral range covered, 
and (2) without it. Thus, the synthetic spectra 
are derived by dividing (1) by (2). Note that 
during the modelling of all spectra, a constant 
wavelength step of AA = 0.05 Awas assumed, 
which ensures errors of the flux in the contin¬ 
uum and in the lines no greater than 0.005% 
of the residual intensity. During the analysis of 
moderate-resolution spectra (FWHM < 5 A), it 
is possible to investigate not only the individ¬ 
ual lines but also broad blends AA > 5 A con¬ 
sisting of lines and bands of many atoms and 
molecules. Therefore, to determine the param- 
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Figure 8. The fit of the isochrone parameters to the observed spectrum of NGC5904 (the profiles of the 
H/3, H 7 , H(5 lines). The three upper panels show the applied model log t = 10.1, Y = 0.30, Z = 0.001 (see 
the corresponding C-M diagram on the left of Fig. 4). The three lower panels show the model log t — 10.15, 
Y = 0.26, Z = 0.001 (the right-hand side of Fig. 4) for comparison. 


eters and abundances of GCs, it is necessary to 
achieve the best fit between theoretical and ob¬ 
served spectra in the whole studied range. As 
the analysis in [22] with S/N > 100 showed, this 
approach allows us to determine the cluster pa¬ 
rameters and the abundances of about 10 chem¬ 
ical elements. Our abundances, determined by 
the strong Ca, Mg, Fe, CH lines and blends, the 
majority of the Na, Al, Ba, Sr lines, and all the 
lines with A > 5300 A, are differential, as the 
empirical oscillator strengths gf, derived in [36] 
by the solar spectrum analysis, are used. For the 
rest of the lines, theoretical values of gf are used 


which can cause, according to the analysis of Shi- 
manskaya et al. [36], a systematic abundance un¬ 
derestimation of up to 0.07 dex. A common mi¬ 
croturbulence velocity ^turb fob all the GC stars 
(Table 3) is derived from the best fit between the 
strong and weak iron lines in the theoretical and 
observed spectra. Figures 10-13 show examples 
of Fe line variations at different ^turb values. Ta¬ 
ble 3 shows the ^turb values derived by us. Note 
that the use of a common value of ^turb fob all the 
cluster stars is not quite proper. However, the 
analysis of the integrated spectra of GCs does not 
allow us to derive a dependence of ^turb on stellar 
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Figure 9. The fit of the isochrone parameters to the observed spectrum of NGC6254 (the profiles of the 
H/3, H 7 , H(5 lines). The three upper panels show the applied model log t = 10.05, Y = 0.30, Z = 0.0004 (see 
the corresponding C-M diagram on the left of Fig. 5). The three lower panels show the model log t = 10.15, 
Y — 0.26, Z — 0.0004 (the right-hand side of Fig. 5) for comparison. 


parameters. Moreover, the main contribution to 
the optical spectra is made by the stars located 
at the MS turnoff point at Tgfj = 6200-7500 K 
and by the brightest asymptotic branch giants. 
For such stars, the microturbulence velocity is 
within the range of ^turb = 1.8-2.4 kms“^. That 
is why the use of a common value of ^turb can¬ 
not cause any significant errors in the abundance 
of chemical elements. The macroturbulence and 
rotation of individual cluster stars are neglected 
in our calculations, as those effects amount to 
less than 30 kms“^ and, consequently, are in¬ 
significant at our resolution. For the same rea¬ 


son, we do not consider the velocity dispersion of 
the stars within the cluster. 

Fitting the continuum of the observed spec¬ 
tra to the theoretical ones is conducted by 
the successive use of the following digital fil¬ 
ters (see the description of MIDAS^ for details): 

(1) determining the maximum spectral intensity 
within a specified wavelength range (5FWHM), 

(2) smoothing filter (running median) with a 
radius of 6FWHM, where FWHM denotes the 
spectral resolution of the spectrograph in A. 

^ http:// WWW .eso.org/sci/software/esomidas/doc/ 
/user/98N0V/volb/ 
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Figure 10. Determination of microturbulence velocity ^turb and Fe, C, Mg, Ca, and Cr abundances for 
NGC6229. The observed spectrum (solid lines) and theoretical spectra (dashed, dashed-and-dotted, and 
dotted lines) calculated for different [X/H] values with constant abundances of other elements (see Sections 

3.2, 4) are shown. 








Relative intensity Relative intensity Relative intensity 


Observ 

t =2.1 km/sec 

nurb 

E =1.8 km/sec 

nurb 


- Observ 

[Fe/H] = -1.9 
-[Fe/H] = -1.8 


vi 11 A 

'':'7 \ ' ■ 


r, < 

-7 ' 


-' mi 'i 


0.90 ^ V 


5120 5160 5200 5240 


5160 5200 


Observ 

[C/H] = -2.05, [0/H] = -2.1 
[C/H] = -1.95, [O/H] = -2.1 
[C/H] = -2.05, [O/H] = -1.9 


- Observ 

[Mg/H] = -2.0 
-[Mg/H] = -1.9 


* 1 . / ' 

0.9 ^ f 




/ \ 

\ f ' 

■ 1 t \ 

1 I 1 

/ '\ 

/ '' 

/ ^ 


4280 4300 4320 4340 


5160 5180 


- Observ 

[Ca/H] = -1.65 
-[Ca/H] = -1.75 

\ 

' / \\ 

V 11 


1/ 

-1//-,-1- 


- Observ 

[Cr/H] = -2.0 
1.00 --[Cr/H] = -1.8 





























19 


Table 3. Comparison of the values taken from the literature [1, 17] and those determined by us using photometric 
and spectroscopic methods: [Fe/H] (in dex), helium abundance Y, age t (in Gyr), a-element abundance 
[a/Fe\ = {[Mg/Fe\ + [Ca/Fe\) /2 (in dex), and average microturbulence velocity for each cluster ^turb in kms“^. 
The parameter HER = [B — R)/{B + V + R) [1] characterizes the number of stars in the blue and red regions of the 

horizontal branch and in the Hertzsprung gap 


Parameter 

NGC 6229 

NGC 6779 

NGC 5904 

NGC 6254 

Ref. 

[Fe/H], dex 

-1.43 

-1.94 

-1.27 

-1.52 

[1] 

HER 

0.24 

0.98 

0.31 

0.98 

[1] 

[a/Fe] 

- 

- 

0.21 

0.26 

[37] 

[a/Fe] 

- 

- 

0.31 

0.39 

[31] 

[Fe/H] 

- 

-2.00 

-1.33 

-1.57 

[17] 

t 

- 

12.75 ±0.50 

11.50 ±0.25 

11.75 ±0.38 

[17] 

Y 

- 

0.25 

0.25 

0.25 

[17] 


-1.74 

-2.35 

-1.33 

-1.74 

ours 

^isochr.+spectra 

12.6 

12.6 

12.6 

11.2 

ours 

^sochr.+spectra 

0.26 

0.23 

0.30 

0.30 

ours 

[F spectra 

-1.65 

-1.9 

-1.6 

-1.75 

ours 

\otl F e] spectra 

0.28 

0.08 

0.35 

0.025 

ours 

^turb 

2.1 

2.1 

2.1 

1.8 

ours 


The contribution of the stars of different types 
to the integrated GC flux is described by the 
example of NGC6229 in Fig. 14. This cluster 
has an HB of an intermediate type, although its 
blue and red ends are densely populated with 
stars. Hot blue HB stars determine the width 
and intensity of the hydrogen lines. These ob¬ 
jects contribute up to 40% of the energy distribu¬ 
tion in the blue region of the investigated spectral 
range. Red giants, dominating in its red region, 
are the main source of molecular bands and lines 
of neutral atoms of heavy elements. Note that 
the molecular lines in the sum spectrum of the 
stars in the red giant branch are much more no¬ 
ticeable than in the spectrum of the MS and sub¬ 


giants branch stars. As a result, the integrated 
spectrum of the whole cluster noticeably differs 
from the spectra of all its components. 

3.3. Abundance Determination 

Variations of [Fe/H], log g, and Tgff allow us 
to calculate the spectra of various cluster stars. 
Summing up these spectra taking into account 
the stellar mass functions produces the total syn¬ 
thetic spectrum of a GG. 

The average chemical abundances in the at¬ 
mospheres of GG stars are determined by match¬ 
ing the observed and theoretical line profiles, line 
blends, and molecular bands. As explained in the 
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Table 4. Derived abundances [X/Fe] (in dex) and their dispersion a (in dex) for NGC6229 and NGC5904. The 

corresponding values from [31] are listed in the last column. 



NGC 6229 

NGC 5904 

NGC 5904[31] 

Elem. 

[X/Fe] 

a 

[X/Fe] 

a 

[X/Fe] 

a 

C 

-0.05 

0.15 

-0.20 

0.10 

-0.48 

0.26 

N 

0.35 

0.25 

0.20 

0.20 

0.68 

0.59 

0 

0.35 

0.30 

0.20 

0.30 

0.15 

0.27 

Na 

0.25 

0.20 

0.20 

0.15 

0.19 

0.26 

Mg 

0.15 

0.15 

0.40 

0.15 

0.33 

0.10 

Ca 

0.40 

0.15 

0.20 

0.10 

0.28 

0.11 

Ti 

0.40 

0.25 

0.20 

0.20 

0.22 

0.10 

Cr 

0.15 

0.20 

0.10 

0.15 

-0.08 

0.19 


Table 5. Derived abundances [X/Fe] (in dex) and their dispersion a (in dex) for NGC6779 and NGC6254. The 

corresponding values from [39] are listed in the last column. 



NGC 6779 

NGC 6254 

NGC 6254[39] 

Elem. 

[X/Fe] 

a 

[X/Fe] 

a 

[X/Fe] 

a 

C 

-0.15 

0.15 

-0.15 

0.10 

-0.77 

0.37 

N 

-0.10 

0.25 

0.25 

0.20 

1.01 

0.45 

0 

-0.20 

0.30 

0.40 

0.30 

0.23 

0.24 

Na 

-0.20 

0.20 

0.25 

0.15 

0.17 

0.27 

Mg 

-0.10 

0.15 

0.05 

0.15 

0.44 

0.13 

Ca 

0.25 

0.15 

0.00 

0.10 

0.33 

0.11 

Ti 

0.10 

0.25 

0.45 

0.20 

0.26 

0.12 

Cr 

-0.10 

0.20 

-0.09 

0.15 

0.001 

0.15 


previous paragraph, the resolution of our spectra 
does not allow us to determine chemical abun¬ 
dances from individual lines. In view of that, we 
adopted such abundances for each element for 
which all the blends (including the lines of this 
element in the studied wavelength interval) are 
best described by the theoretical spectrum. 


of the elements whose lines and/or molecular 
bands dominate in the spectra which allows 
us to obtain [X/H] values with an accuracy of 
A[X/H] ~ 0.1-0.2 dex. These elements are Fe, 
Ca, Mg, C. The line profiles of Ti, Cr, Co, Mn, N 
are analyzed by varying their abundances using 
previously derived fixed abundances of the el¬ 
ements from the first group. The accuracy of 


First, we determine the chemical abundance 
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Figure 14. The contribution of different types of stars to the integrated GC radiation. The example of 

NGC6229 is considered. 


chemical abundance estimation for the elements 
of the second group, the lines of which are much 
weaker and more blended than the lines of the el¬ 
ements from the first group, is about 0.2-0.3 dex. 
When we use moderate-resolution spectra, the 
individual elements do not have any noticeable 
details in the observed wavelength range, but 
they influence the molecular and ionization equi¬ 
librium of other elements. Oxygen belongs to 
such elements; as its abundance increases, a part 
of atomic carbon participates in the formation of 
the CO molecule which reduces the intensity of 
the CN and CH molecular bands. The elements 


Al, Si, V, Ni influence the electronic equilibrium 
in stellar atmospheres only slightly, and also do 
not have observable lines at moderate resolution. 
Na should also be included in this group, as its 
resonance lines are heavily distorted by interstel¬ 
lar lines and are unsuitable for analysis. 

4. RESULTS 

This is the first time that the chemical abun¬ 
dances for the clusters NGC 6229 and NGC 6779 
have been derived. Tables 3-5 and Figs. 10-13 
show the final results for them and for the com¬ 
parison clusters NGC 5904 and NGC 6254. Ta- 
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ble 3 summarizes the GC fundamental param¬ 
eters: age, helium abundance, iron [Fe/H] and 
a-element abundance averaged over the GC. Ta¬ 
bles 4 and 5 show the abundances of individual 
chemical elements. Figs. 6-9 show the results 
of the comparison of the observed and theoret¬ 
ical H/3, H 7 , H5 line profiles, calculated based 
on the isochrones with the parameters derived 
from the photometric data (see above). Figs. 10- 
13 show the profiles of the lines and molecular 
bands of different chemical elements that were 
used to estimate their abundances. For each ele¬ 
ment, two-three variants of the line profiles with 
different abundances are presented for reliability 
assessment and accuracy check. It is necessary to 
note that when fitting the theoretical spectrum 
to the observed one, we assume that the light 
from GC stars in the main evolutionary stages is 
adequately represented in these two spectra. It is 
also assumed that the observed spectrum is not 
distorted by background objects. In reality, the 
last condition is not always met, especially in the 
case of GCs close to the Galactic plane. Obser¬ 
vations with different slit orientations allow us to 
eliminate the distortions caused by background 
stars in most cases. The spectra of test GCs 
from Schiavon et al. [3] were derived using the 
long-slit scanning method, and, thus, they may 
contain such distortions. 

As was mentioned in the Introduction, age, 
[Fe/H], and V were estimated using deep pho¬ 
tometry results from the literature. Table 3 
shows the summary of those estimates. Gener¬ 


ally, our results of the above-mentioned estima¬ 
tion of these parameters, derived by means of 
cluster spectrum analysis, agree with the ones 
from the literature. However, there are a num¬ 
ber of discrepancies. Our age estimate for the 
well studied GC NGC 5904 exceeds the value de¬ 
rived in [17] by 0.9 Gyr, and the helium abun¬ 
dance is excessive. These parameters are de¬ 
termined by the HI line profile analysis, as ex¬ 
plained in Section 3.2. Our [Fe/H] values, de¬ 
rived using the spectroscopic method, are in most 
cases smaller by 0 . 2 - 0 .3 dex than those from the 
literature, derived by the C-M diagram analy¬ 
sis. It is unlikely that these discrepancies are 
caused by background stars caught in the slit. 
The projected vicinities of the investigated clus¬ 
ters mainly contain stars of the Galactic disk 
with close-to-solar metallicities. We have not 
detected any emissions in the observed cluster 
spectra. Thus, the noted differences are likely 
associated with the distinctions of the applied 
methods and models. When analyzing the pho¬ 
tometric data, the errors of color indices and dis¬ 
tances to the clusters, and also the use of theo¬ 
retical isochrones from different authors produce 
random and systematic shifts in the derived pa¬ 
rameters [38]. The use of different atomic and 
molecular data sets, different model atmospheres 
for the stars, and different principles of modelling 
simple stellar populations leads to differences in 
the abundance estimates and cluster evolution¬ 
ary parameters. 

Tables 4 and 5, and Figs. 15 and 16 show the 
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comparison of our abundance estimation results 
for different elements and the results from [39] 
and [31] for the well studied clusters NGC5904 
and NGC6254. The abundance differences are 
greater than the given errors only for G and N. 
A possible explanation is that we analyze the av¬ 
erage chemical abundance of the whole cluster 
and not only of the high-luminosity stars, as the 
majority of the authors do. The products of the 
GNO cycle could be brought up to the surface 
of these stars, and this may cause an increase in 
the nitrogen abundance and a decrease of that of 
carbon. Note that the considerable contribution 
of red giants to the cluster spectrum probably in¬ 
creases the distortions of the estimated C and N 
abundances, but their amplitude should be much 
smaller than in the investigation of individual 
stars. There are also significant differences in the 
Mg and Ca abundances for NGG 6254. Gonsider- 
able variations of C, N, O, Mg, Na, and A1 from 
star to star were found for a number of massive 
Galactic GCs [37]. In one and the same clus¬ 
ter, one can often observe two or more stellar 
populations with different chemical abundances 
forming different sequences in the C-M diagram. 
This is explained by the mixing of matter after 
the GNO cycle. The Mg abundance variations 
after such processes are expected to be around 
AMg < 0.3 dex [39]. Note once again that the 
estimates from the literature are based mainly on 
the investigations of red giants, which are not the 
true indicators of GC initial abundances because 
of the mixing of matter, to which their atmo¬ 


spheres are liable (the first dredge-up). There¬ 
fore, more detailed theoretical and observational 
investigations are necessary in order to check how 
much the cluster abundance variations depend on 
the physical characteristics of the studied stars. 

Almost all our estimated abundances agree 
with the theoretical calculations of the models 
of the chemical evolution of the Galaxy and also 
with the average values for dwarfs at the given 
metallicity (see [40-42]). However, the scatter of 
the observed values, with which model character¬ 
istics are compared, is often large. 

Carbon is produced in the triple a-process 
with hydrostatic He burning in stars of various 
masses. The average abundance of C in halo and 
disk dwarfs varies little with time (does not de¬ 
pend on [Fe/H]): [C/Fe] ~ 0 [41]. We have de¬ 
rived similar values. 

The nitrogen abundance decreases with 
metallicity as massive stars do not produce pri¬ 
mary N. With the [Fe/H] increase, the produc¬ 
tion of the so-called secondary N by moderate- 
mass stars increases until the Fe production 
by la supernovae compensates this process at 
[Fe/H] ~ —1. The theoretical [N/Fe] values 
vary from 0 to —0.5 in the following range: —2 < 
[Fe/H] < 1. Our estimates are closer to the the¬ 
oretical ones than the data from the literature 
for the giants in NGG 5904 and NGG 6254. 

Oxygen, as the majority of other a-elements, 
is produced by massive stars only. [O/Fe] ~ 0.5 
at [Fe/H] < 1. Our values correspond quite 
well to those calculations, except for NGG 6779, 
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A comparison of the chemical abundances in NGC 6229 and NGC 5904 according to the results 
of this paper. The open circles denote the data from [31]. 



Figure 16. A comparison of the chemical abundances in NGG 6779 and NGG 6254 according to the results 
of this paper. The open circles denote the data from [31]. 
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for which we obtained lower abundance esti¬ 
mates for the other a-elements as well. There 
were no emissions or noticeable distortions from 
background stars detected in the spectra. How¬ 
ever, one cannot exclude the latter possibility for 
NGC 6779, which is closer to the plane of the 
Milky Way than the other three clusters. 

The evolution of the abundance of magnesium 
as a function of metallicity is impossible to ex¬ 
plain only by its formation in massive stars ([41] 
and references therein). The theoretical [Mg/Fe] 
values in the range of — 2 < [Fe/H] < 1 are usu¬ 
ally a little lower than the observed ones. Our 
results agree better with the theoretical values. 
However, the scatter of the observed data is large. 
On the whole, one can conclude that the [Mg/Fe] 
values for the selected clusters except NGC 6779 
are within the normal range. 

The value of [Ca/Fe] is in close agreement 
with theory for all the selected clusters except 
NGC 6254. We estimated the abundance of this 
element not only from the Cal A 4226 A reso¬ 
nance line but also from the H and K Call lines, 
the profiles of which vary greatly depending on 
the selected isochrone, as these profiles are de¬ 
termined not only by the Ca content but also by 
age, y, and [Fe/H]. Thus, an error in the [Ca/Fe] 
estimate is possible, but much less probable than 
the possible errors for the elements heavier than 
Ca, which have no broad spectral features. 

Tables 4 and 5 and Figs. 15 and 16 show 
that, on the whole, the chemical abundances for 
NGC 5904 and NGC 6254 are similar to the ones 


for NGC 6229 and NGC 6779. The published val¬ 
ues of [Fe/H] for NGC 6779 and NGC 6254 are 
smaller than the same values for NGC 6229 and 
NGC 5904. This fact, in view of the similar¬ 
ity of the horizontal branches of NGC 6229 and 
NGC 5904, and also of NGC 6779 and NGC 6254, 
has been a reason for associating the clusters of 
higher metallicity with the young Galactic halo, 
and clusters of lower metallicity—with the old 
one [2]. Following our results (see Table 3), there 
are no two objects out of the four investigated 
that are totally similar in all parameters. The 
similarity in the chemical abundances is indica¬ 
tive of the fact that they belong to the old Galac¬ 
tic halo, in which nucleosynthesis proceeded ac¬ 
cording to the main evolutionary stages under 
the influence of SNH and SNIa supernovae. In 
order to associate the clusters with one or an¬ 
other Galactic subsystem, a more detailed anal¬ 
ysis of the chemical, kinematic, and structural 
characteristics is necessary. 

5. CONCLUSIONS AND SUMMARY 

The method we used for modelling and anal¬ 
ysis of the integrated spectra of the clusters was 
previously discussed in [22]. In this paper we ex¬ 
panded it by including the investigation of the 
cluster C-M diagrams and their detailed com¬ 
parison with the theoretical isochrones. Now, 
as a result, this method combines the possibil¬ 
ities of using both the spectroscopic, integrated 
along the slit data, and photometric data on the 
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stellar populations of the clusters. We use the 
results of deep stellar photometry from the lit¬ 
erature to select an isochrone that best fits the 
cluster C-M diagram. The summing-up of the 
stellar synthetic blanketed spectra, calculated us¬ 
ing model atmospheres, is conducted using the 
stellar luminosity function of Chabrier [33]. The 
matching of the theoretical and observed clus¬ 
ter spectra is carried out iteratively. Initially, 
the distribution of cluster stars by mass, radius, 
and log g is preset using the literature data on 
the theoretical isochrone which fits best the ob¬ 
served C-M diagram. For the specification of 
stellar ages and helium abundances, the profiles 
of the hydrogen lines are analyzed, because they 
practically do not depend on other parameters. 
The abundance of iron is then varied, the lines of 
which prevail in the optical spectrum even at a 
low metallicity. Upon achieving the closest agree¬ 
ment in these three parameters, the abundances 
of the a-elements Mg, Ca, and C are derived. 
The molecular bands and lines of these elements 
are the most significant in the spectrum. For 
the determination of the abundances of other el¬ 
ements, complex blends consisting of many lines 
are used. As a result, the abundances of about 
10 elements are estimated. Note that the use 
of a considerable number of lines with empirical 
gf values allows us to obtain differential abun¬ 
dances (i.e., bound to the gf values of solar spec¬ 
tral lines) of Ca, Mg, Fe, C, and, in some cases, 
other elements. 

In this paper, we determined the ages. 


specific helium abundances Y, and abundances 
of Fe, C, N, O, Na, Mg, Ca, Ti, and Cr for the 
following four Galactic GCs using the developed 
method: NGC6229, NGC6779, NGG5904, 

NGC 6254. The chemical abundances for the 
first two objects were estimated with a spec¬ 
troscopic method for the first time. According 
to our estimates, these GGs turned out to 
be about 1 Gyr younger than NGC 5904 and 
NGC 6254, which are accepted in the literature 
as their analogues. The helium abundance 
that we determined for the last two clusters is 
higher (Y = 0.30) than the one in the literature. 
The clusters NGC 6229 and NGC 6779 show 
the values Y = 0.23 and Y = 0.26 respec¬ 
tively, which are closer to the value of the 
primordial helium abundance in the Universe 
Y = 0.25 [43]. The a-element abundances 
in the metal-poor {[Fe/H] ~ —1.6) NGC6229 
and NGC 5904 turned out to be high 
{[a/Fe] = 0.28-0.35) as in the old massive 
globular clusters of the Galaxy; massive type II 
supernovae contributed greatly to their chemical 
enrichment. However, the a-element abundances 
in the investigated clusters with low central 
density, NGG6779 and NGC 6254, turned out 
to be quite low {[a/Fe] = 0.08,0.025), while 
their metallicity was low {[Fe/H] = —1.9, —1.7). 
Such [a/Fe\ values are not typical for low 
metallicity Galactic clusters. For this reason, 
the search for a dependence of the chemical 
abundance of GGs on their structure and 
kinematic characteristics can be the subject 
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of future investigations. A comparison of the 
derived abundances and the literature values for 
the well studied GCs NGC5904 and NGC6254 
showed that our measurements are consistent 
with the corresponding average estimates from 
the literature for cluster giants with the average 
accuracy of about 0.15 dex, except for C and N. 
This is probably due to the difference between 
the average chemical composition in the stellar 
atmospheres for the whole cluster, derived in 
this paper, and the chemical composition of 
the atmospheres of high-luminosity stars listed 
in the literature. The products of the GNO 
cycle can be transported to the surface of such 
stars, and their study does not provide accurate 
information on the primary chemical abundance 
of GGs. 

The values of [Fe/H], Y, and age which we de¬ 
rived by spectroscopic analysis allow us to eval¬ 
uate the true correlation between the observed 
C-M diagram and the corresponding theoretical 
isochrone. At present, there is no possibility of 
combining the deep stellar photometry data with 
the spectroscopic methods of estimating the el¬ 


ement abundances and evolutionary parameters 
when analyzing the characteristics of stellar pop¬ 
ulations in extragalactic GCs. Thus, our investi¬ 
gation provides essential data for understanding 
the characteristics of GCs in remote galaxies. 
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